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ABSTRACT 



Context. Properties of candidate stars, forming out of molecular clouds, depend on the ambient conditions of the parent cloud. We 
present a series of 2D and 3D simulations of fragmentation of molecular clouds in starburst regions as well as clouds under conditions 
in dwarf galaxies, leading to the formation of protostellar cores. 

Aims. We explore in particular the metallicity dependence of molecular cloud fragmentation and the possible variations in the dense 
core mass function, as the expression of a multi-phase ISM, due to dynamic and thermodynamic effects in starburst and metal-poor 
environments. 

Methods. To study the level of fragmentation during the collapse, the adaptive mesh refinement code FLASH is used. With this code, 
including self-gravity, thermal balance, turbulence and shocks, collapse is simulated with four different metallicity dependent cooling 
functions. Turbulent and rotational energies are considered as well. During the simulations, number densities of 10 8 -10 9 cm -3 are 
reached. The influences of dust and cosmic ray heating are investigated and a comparison to isothermal cases is made. 
Results. The presented results indicate that fragmentation increases with metallicity, while cosmic ray and gas-grain collisional 
heating counteract this. We also find that modest rotation and turbulence can impact the cloud evolution as far as fragmentation is 
concerned. In this light, we conclude that radiative feedback, in starburst regions, will inhibit fragmentation, while low-metallicity 
dwarfs also should enjoy a star formation mode in which fragmentation is suppressed. 
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1. Introduction 

The initial mass function (IMF) in our local neighborhood 
^ | . is observed to be a power-l aw function n icely following a 
1 Salpeter slope. Recent studies (lFigedl2Q05blfal : ISabbi et al.ll2007l: 
\ lElmegreen et al.ll2008h confirm this universality at low and high 
r* ■ masses. When we turn to observations of more radical environ- 
' ments like the Galaxy center, a 'universal' shape is again con- 
K> „ firmed and only hints f or different IMFs appear. Measurements 
;_j of abundance patterns dBallero et al.ll2007l: ICunha et al.ll2008l) 
indicate a flatter IMF. O bservations of the Arches cluster init ially 
showed a shallow IMF dFiger et al.lll999l: IStolte et al.ll2002l) and 
a turn-o ver mass, and thus a typical m ass, at a few solar masses 
(6-7M^, IStolte et al]l2005h . However. [Kim et al. I (I2006I) do de- 
tect low mass stars with their deep photometry and merely find 
the presence of a local bump at ~6.3M Q and a s omewhat shal- 
lower IMF (r = -1.0 to -l.D. lDib etall (12007b confirm these 
results and claim a top-heavy IMF. Again, there appear to be per- 
fectly reasonable expla nations for the higher mass st ars detected 
in the Arches cluster as lPortegies Zwart et al.l (120071) argue with 
their idea of a still collapsing cluster core. Of course, the idea 
that the conditions and the environment play an important role in 
shaping the IMF is worth studying. Lookin g at the mass-to-light 
ratios of ultra-compact dwarf galaxies, iDabringhausen et all 
( 2009) conclude from their models that it is most likely due to 
a top-heavy stellar IMF. If so, then they attribute this to feed- 
back effects induced by massive stars. Still, these studies do 



not benefit from resolved stellar population data. One can asso- 
ciate the conditions and the environment o f dwarf galaxies with 
those in the early Universe dHunted [20081: ISalvadori & Ferraral 
2009). It is also believed that the primordial, high redshift 



and zero metallici ty, IMF must hav e been top-heavy, ~100Mq 
1 [20001: lOmukai & Pallal I2001L iBromm etall 120021: 



(lAbel 



et a 



Abe l et al 



2002; Omukai&Palla 2003; 



Hence, it appears that the IMF is universal, but that it is worth- 
while to explore strong environmental (feedback) influences, 
e.g., starburst and dw arf galaxies, to understand the IMF better 
dKlessen et al.ll2007l) . Since star formatio n, fragmentation and 
the IMF all depe nd on initial conditions (ISkillman et al.l l2003: 
iGreif et al.l l2008). the question is how strong the impact of these 
initial conditions is. Also, how different should the environment 
be, or how strong the feedback, to cause a significant deviation 
from a Salpeter IMF? 

Molecular clouds are the birth grounds of stars. These 
clouds are formed (rapidly) when interstellar gas and dust 
gather and cool down under the i nfluence of hydrodynam- 
ics, ther modynamics an d gravit y dBallesteros-Paredes et al.l 
19991: iHartmann et aD I2OOII: iKlessen et alJ 

Dobbs & Bonnelll 120081: 



.Har tmann e t al.. 
Heitsch & Hartmannl 120081: 



2008)). Molecules form in the gas phase and on dust grains , 
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and allow low temperatures to be reached (lOmukai et alj|2 005 ; 
ICazaux et al.ll2005l: ICazaux & Spaansll2009l: iDulieu et al.ll2009h . 
Dense regions, like the center of the cloud, are able to form 
many stars , with a modest efficiency (-5%: lElmegreen| [2000: 
IClark et al.l l2005h . The distribution of the dense regions is 
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thought to b e the precursor of a stellar initial mass function 
dMotte et alJl!998t Motte & Andrdl2Q0H iGoodwin et alJl2008h 
ISimpson et al.l 120081) . and perhaps the stellar IMF is directly 
linked to this core mass function as it seems for the Pipe nebula 
(lLada et al.ll2008l) . Following a typical molecular cloud from its 
infancy to a mature state, for different starting conditions, gives 
insight into the processes that determine stellar masses and is 
the focus of this work. Specifically, we consider low metallicity 
environments (dwarfs) and warm dusty systems (starbursts). 

In the next section we present the initial conditions used 
for our 38 simulations and present the results in the following 
section. Twenty of these simulations form the basis of this re- 
search. They comprise 4 different metallicities with 5 combina- 
tions of turbulent and rotational energies each. We add to these 
basis simulations comparison studies that include extra heating 
sources and cases where the primary heating and cooling terms 
are removed. 

2. Simulations 

2.1. Numerical Model and Simulation Setup 

In order to model a collapsing molecular cloud, we solve the 
hydrodynamical equati ons using the adapti ve mesh refinement 
(AMR) code FLASH dFrvxell et aD l2QQQ|). FLASH uses the 
PARAMESH library (lOlson et alJll999b iMacNeice et alJl2000h 
for grid refinement and parallelization. Hydrodynamic equa- 
tions are solved usin g the p iecewise parabolic method (PPM, 
IColella & Woodwardl dl984b ). 

For our 2D simulations we use a grid with a maximum res- 
olution that amounts to 4096 2 cells when completely refined, 
which is a refinement level of 10 in FLASH terms. To mini- 
mize computational demands, we utilize the power of adaptive 
meshes by initializing the simulation at refinement level 5, i.e., 
128 2 cells, and refine with a density thresh old ofn = 600 cm" 3 i n 
order not to violate the Truelove criteria dTruelove et al.ll 19971) . 
This density was found to be the optimal balance between fol- 
lowing the physics of cloud collapse and computational demand. 
The most dynamic 'central' area, however, always starts with 
the highest refinement to maximize precision and minimize res- 
olution dependent effects in this region. The border resolution 
between the refinement levels is set-up in such a way that it is 
not a sharp transition between the minimum and the maximum 
refinement. 

We are aw are of the problems th at might arise if resolution 
is not sufficient. iTruelove et al.l (1 19971) state that the Jeans length 
should at least be resolved by 4 cells, i.e., N/ > 4 in order to 
ensure that the collapse is of a physical rather than numerical 
nature. Another pitfall is that the transfer of gravitational energy 
to rotational energy can lead to unphysical values if the initial 
resolution is poor, i.e., t he number of initial cells N; < 2145 
(iCommercon et al 3 1200% . We are well beyond these limits. Our 
minimum number of cells, extending beyond the central region, 
is at least N; > 16384 cells. The Jeans length is resolved by 
more than N 7 > 77 cells. For even greater precision, we keep 
the tree structure of the top level hierarchy intact in FLASH. 
In other words, we initiate all simulations with one top-level 
block. Every simulation is continued up to twelve dynamical 
timescales, which is computed as 

fdyn = ^3n/32Gp. (1) 

Where G is the gravitational constant and p is the mass 
density. 



We do not include magnetohydrodynam ic (MHD) effects 
in our models. iTillev & Pudritzl (120051 120071) argue that these 
should not matter for the core formation as much as the turbu- 
lence does. However, MHD certainly is an important aspect, see 
iBanerjee et al.l (2007, 2008), that we will investigate further in 
the future. 

2.2. initial Conditions 

We start the basis simulations with four metallicities Z/Z© = 
1, 10" 1 , 10" 2 , 10" 3 with respect to solar metallicity and take Z = 
Z© as a benchmark value. The dependence of the model out- 
comes on metallicity is the focus of this work. 

For all simulations, we impose Gaussian turbulence with a 
charact eristic dispersion of cr = 2 km/s, compara ble to observed 
values (Fal garone et al]l2001c ICaselli et lri]l2002l) and of the or- 
der of 0.3 ^GM/R. We do not assume the clouds to be rotation- 
ally supported but rather that turbulent motions, following the 
Larson laws, dominate. We adopt a Larson like sc ale relation fo r 
the turbulence of the form cr = 2 (R/lpc) 5 km/s (lLarsonll98l|) . 
iBonnell et al.l (120061) also show that molecular clouds exhibit a 
cr oc R 5 velocity dispersion law. The turbulence is introduced as 
part of the initial conditions of the cloud and is not maintained by 
external forcing. The virialization of the gas that falls into local 
potential wells does maintain a level of super-thermal random 
motions for several dynamical timescales. 

We run each metallicity condition with five rotational en- 
ergies, i.e., /3 = 10 _1 , 10" 2 , 10" 3 , 10" 4 , 0, where fa is the initial 
ratio of rotational to gravitational energy. A Keplerian rotation is 
initiated around the z-axis multiplied by the ratio fio. We define 

I co 1 v 2 r 

* = ~&Wr = ~GM (2) 

Here, / is the moment of inertia, co is the angular velocity, 
v is the radial velocity and M and r are the total cloud mass and 
its corresponding radius. Studies show that a typ ical ratio for 
molec ular cloud cores is of the order fio = 0.02 (ICaselli et al.l 
120021) . These authors have observed 57 cloud cores and find en- 
ergies within the range of fio « 10" 1 - 10" 4 . Cosmological simu- 
lations indicate that star form ing host clouds have rotation ener- 
gies o f the order of fa - 0.1 (iBromm et al.ll2002t lYoshida et all 
l2006h . 

We simulate an initially spherical cloud with a radius of 10 
pc within a box that is periodic in both gravity and spatial bound- 
aries. Given the larger box size of 26 pc and due to symmetric 
initial conditions, we do not expect nor see that the periodicity 
is an issue. We initializ e the cloud with a mean constant den- 
sity of n 10 3 cm -3 (lFrieswijkll2008l) . since it is not likely 
that a well developed density profile exists before virialization. 
These initial conditions mean that we start with an initial cloud 
mass of 28.3xl0 3 solar masses. The initial gas temperature is 
put at 100 K throughout the cloud, s imulating starburst regions 
(ISpaans & Sii3l2000b lOtt et al.ll2005l) . 

2.3. The Simulations 

We present twenty simulations, with flat initial density profiles, 
where we focus on the effect of metallicit y on cloud fragm en- 
tation. In a similar fa shion as the work by iMachidal d2008) and 
iMachida et al.l (120091) . we also investigate the importance of 
rotation and turbulence in these. We present 15 more simulations 
with the same setup to test and compare specific physical cases. 
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Metallicity: Four sets of five simulations are done in order to 
follow the collapse and fragmentation of metallicity dependent 
cooling functions for turbulent, rotating and non-rotating cases. 
These simulations have an initial cloud temperature of 100 
K, which represents the high temperatures within the active 
environments where the molecular clouds are formed. 
Dust and CRs in starbursts: Five solar metallicity runs with 
additional heating sources; cosmic rays and dust. In massive 
star forming regions prominent in starbursts there can be a lot of 
warm dust affecting the thermal balance of the cloud. Dust can 
act as a cooling mechanism when the gas temperature is higher 
than the dust temperature. Once the temperature of the gas is 
lower than the dust, and this is the case for photon dominated 
regions (PDRs) with Ay > 5 mag, it will become a source of 
heating. 

Cosmic rays can ionize and heat the gas, but their contribution 
is generally thought to be low. In starbursts, the impact of 
CR-heating is stronger due to the increased production by 
massive stars, i.e., supernovae. 

Isothermal at 25 K: Five runs kept the temperature nearly 
constant at a starting cloud temperature of 25 K. Representing 
an environment where the temperature is higher than what we 
observe in the Milky Way. For these simulations, the temper- 
ature is not allowed to go below 25 K and the compressional 
heating term is removed. The results are compared to cases 
where the thermal balance is incorporated. 
Isothermal at 10 K: A second series of near-isothermal runs 
are initiated at T=10 K in order to compare to Milky Way 
conditions and to check the impact of the isothermal condition. 



2.4. Chemistry, Cooling and Heating 

In our simul ations we use a de t ailed cooling function cre- 
ated with the iMeijerink & Spaansl (120051) code. The metallicity 
dependent cooling includes fine-structure emissions from car- 
bon and oxygen as well as molecular lines from species like 
CO and H2O. The level populations have been corrected for 
LTE effects above 10 3 cm -3 . The chemistry i ncludes gas phase 
and grain su rface formation of H2 and HD (iCazaux & Spaansl 
l2004ll2009b and line trapping due to the initial column (N H ~ 
1 22 - 5 cm" 2 and cr = 2 km/s ) through the cloud, as described 
in lMeijerink~& Spaans ! (120051) and using t he mu lti-zone escape 
probability method of iPoelman & Spaansl (120051) . Line trapping 
during cloud collapse is ignored (see the discussion). A back- 
ground radiation field has bee n included that conforms to 0. 1 Go 
(in units of the Habing field. lHabindll968) in the dwarfs and 
30 Gp for the starbursts (ISpaans & Normanl 19971: iKlessen et al.l 
12007b . see figureffl 

At large column densities, cosmic ray heating can play an 
importa nt role. For our heating by c osmic rays, we use the defi- 
nition of IMeijerink & Spaansl (120051) as presented in equation [3j 
which scales linearly with the density. 



T C R = 1.5 x 10" 11 jr n(H 2 ) erg cm" 3 s" 



(3) 



We have taken the cosmic ray ionization rate per H2 
molecule f as 9-10" 16 s" 1 . This corresponds to a starburst en- 
vironment with a background star formation rate of the order of 
30 M Q per year, with the appropriate scaling for dwarf galaxies 
(ISpaans & Silkll2005b . 

Besides the heating by cosmic rays, heat can be transferred 
through collisions between dust and gas. The gas-dust tempera- 
ture difference is important here. We add gas-grain collisional 




2.5 3.0 
logT(K) 

Fig. 1. Four cooling curves by Meijerink & Spaans correspond- 
ing to four metallicities Z = Z Q , 10 _1 Z o , 10" 2 Z o and 10" 3 Z o . The 
cooling rates (A cr ) correspond to a number density of 10 4 cm -3 . 



heatin g to o ur s imulations as given by iHollenbach & McKed 
dl989Lll979b and lMeiierink & Spaansl (12005b : 



r co n = 1.2 xio -3 V 



1000/ \ a 



100 A 



(4) 



, -75 

5| tT 



(Td - ?k) erg cm 3 s l . 



In this equation, Tk is the gas kinetic temperature, n is 
the average number density and is the minimum grain 
size that is taken as 10 A We k eep the grain size distribution 



fixed (MRN . iMathis et afll 1977b. The dust temperature we use 

"TUT 



Th = 39 K (fWiklind & He nkel||200lb, is for a starburst region 
(IMeijerink & Spaansll2005b and is kept constant. The heat trans- 
fer between gas and grains will act as a cooling mechanism for 
the gas when it has a higher temperature than the dust. 

Because line trapping beyond a column density 10 22 5 cm -2 
is not taken into account, the temperatures presented in this pa- 
per are lower limits. Above ~10 45 (Z o /Z) particles per cm -3 , 
heating by dust dominates. In this case, line trapping becomes 
less important because the n 2 dependence of the gas-dust cou- 
pling prevails over the n dependence, or weaker, of line emission 
in LTE. The fine intricacies of cooling and heating together with 
gravity determine the main results of this paper. 



3. Results 

First, we perform three 3D simulations at refinement levels 9 and 
10 (2048 3 and 4096 3 cells). We do these simulations to see if the 
extra dimension produces significant differences in the results as 
compared to 2D. Due to one less degree of freedom, 2D simula- 
tions can result in a higher merger rate and thus fewer fragments 
with higher masses in the end. The three 3D turbulent simula- 
tions were carried out for modest rotational energies, = 10" 1 , 
fio = 10" 3 and one non-rotating initial condition, to avoid exces- 
sive influence of rotational motions on the outcomes. These runs 
are slightly lower in resolution compared to the 2D simulations 
and we follow them up to 4 t^ yn . In the two rotating cases, we see 
a disk forming within a time of the order of a free fall timescale 
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t « £ dyn , see figure [2J After r dyn , one effectively has a 2D sys- 
tem. The non-rotating simulation, however, does not exhibit sig- 
nificant flattening. Nevertheless, we do observe that the results, 
even for this case, is very similar to its 2D counterpart for solar 
metallicity. It has been verified that fragmentation starts on the 
same timescales as for the 2D simulations, t > t& yn . The number 
of fragments is somewhat higher while their masses are com- 
parable to their 2D counterparts. We emphasize that metallicity 
effects have our interest and feel that 2D simulations suffice to 
capture any salient features. 

The final stage of every simulation is analyzed for the num- 
ber of fragments (cores) and their masses. The results of all basis 
simulations in metallicity Z and initial turbulent and rotational 
energy, are shown as density images after 12 theoretical dynam- 
ical timescales in figure [3] 

We find the fragments by using an al gorithm called iso - 
volume within the visualization tool Visit ( Child s et al 
A given density threshold value lets the algorithm select the 
connected components above this threshold. Putting a density 
cap at 10 5 cm" 3 t o select the dense cores (lAlves etal] 12001: 
ICaselli et al.ll2008h . we were able to count the number of frag- 
ments and determine their masses. The choice of threshold den- 
sity matters for the number of fragments one derives and their 
masses. However, we find that there is no significant change in 
the number of fragments for a spread of one order of magnitude 
in our threshold density, while the fragment masses do decrease 
with increasing threshold density. Clumps and fragments can 
also be identified from their deep pot ential wells like the bound 
p-cores method by lSmith et al.l (l20Q9h . The fragment sizes range 
between r=6,000 AU and 84,000 AU. A complete overview of 
the fragments and their masses can be found in table [T] 

3.1. Metallicity and Multi-Phase Structure 

We evaluate the temperature inside the cores and that of the sur- 
rounding gas of every simulation. In most cases, the cooling 
function is able to cool the gas efficiently. Even though there 
is compressional heating, the gas can cool down to temperatures 
as low as 10 K. The difference lies in the rate with which each 
model cools down or heats up. For the various metallicities, we 



present temperature plots in figure ffl The images are represen- 
tative for every value of J3q. The density is overplotted with con- 
tours in these images to highlight the high density regions, and 
specifically, to show the location of the dense cores. The Jeans 
length is illustrated on the lower left side of the imaged 

We find that the temperature increases inside the fragments 
as the metallicity decreases. The regions where the density is 
higher than its surroundings have temperatures of the order of 
10 K on average for Z = Z© and 20 K for Z = 10" 3 Z o . We see 
the same trend for the lower density gas, in that the amount of 
cool gas, i.e. gas with temperatures lower than 50 K, decreases 
with Z. For a gas density above 10 4 cm -3 , and a gas temperature 
of more than 100 K, fine- structure cooling of [OI] dominates 
together with mid-J CO emission and water at the 20% level. 
Below 100 K, low- J CO line emission is the main coolant with 
contributions from water on the 10% level. Diffuse gas below 
10 4 cm -3 is mainly cooled by [CII]. 

In figures [5] and [6l we present phase diagrams, gas temper- 
ature versus density, that highlight the response of the gas to 
metallicity driven cooling and additional heating processes. It is 
immediately obvious that a cold dense phase is supported, but 
also that more warm gas exists at lower metallicity or higher 
dust temperatures. In figure [5j one immediately notices that for 
high metallicities, the slope dlog(T)/dlog(n), is steeper between 
1 and 10 3 -10 4 cm" 3 and that the cooling time is shorter. We 
also note that the maximum density achieved is smaller at lower 
metallicities. The diagrams show that a two-phase medium ex- 
ists with preferred densities of around 1 and 10 4 cm~ 3 . This is 
independent of the specific grid structure. 

Above 10 4 -10 5 cm~ 3 , compressional heating starts to dom- 
inate the thermal balance and raises the temperature. The in- 
crease due to adiabatic heating is weakened for higher metal- 
licities, while for lower metallicities the temperature can rise up 
to 1000 K, if collapse is not halted completely. For low turbu- 
lent/rotational energies and in the case of zero rotation, there is 
a distinct curvature at a single density present. This is where the 



1 Note that the diameter of this circle corresponds to twice the Jeans 
length, which we refer to here as the Jeans circle. For an object to be 
able to fragment, it must have a size of at least twice the Jeans length 
and thus, must be larger than the Jeans circle. 




Fig. 3. Density plots of the final stage (t = 12^ yn ) of every simulation of metallicity against rotational energy. Orange to red 
depicts high density, typically higher than 10 4 cm -3 and up to 10 9 cm -3 . Light blue portrays the densities smaller than 1 cm -3 . 
Top to bottom: Decreasing in metallicity, Z = Z Q , 10 _1 Z o , 10" 2 Z o , 10" 3 Z o . Left to right: Decreasing in rotational energy J3q = 
10" 1 , 10" 2 , 10" 3 , 10" 4 , 0. All images have boxsizes of 4x4 parsec. 



temperature rises quickly due to compression so that the inef- 
ficient cooling cannot compensate for this fast increase. At this 
point, the object is thermally supported until it reaches higher 
temperatures where the cooling gets increasingly more effective. 

There are always accretion shocks with high temperatures in 
their wake. This is important for the clouds evolution and frag- 
mentation. We find shocks in low density gas and occasionally 
near the fragments. They consist of thin layers of slightly over- 
dense gas with high temperature and a steep velocity gradient. 
Shock waves are more prominent when the collapse is rapid. 
Shock heating can dominate under such conditions. 



3.2. The Effect of Metallicity, Turbulence and Rotation on 
Fragmentation 

We notice the immediate result that higher metallicity and 
higher rotational energy leads to more fragmentation, also found 



bv lClarketaD (12008^ sim i lar in- dept h work on this has als o 
been performed bv iMachidal (120081) and iMachida etal.1 (12009b . 
at higher densities. 



From a metallicity perspective, a clear trend is visible that 
at higher metallicities, there are more fragments and the average 
masses of the fragments are generally lower. The last statement 
is illustrated in figure [71 The total fragmented mass, if we sum 
all the fragment masses, increases with the amount of fragmen- 
tation, thus with metallicity. This is caused by the fact that the 
fragments accrete mass over time and more fragments accrete 
more mass. The compressibility of the gas also depends on the 
metallicity. One may understand this using a polytropic equation 
of state of the form P oc p?, for an ideal gas P oc pT. Increased 



6 S. Hocuk & M. Spaans: Fragmenting Clouds 

Table 1. Analysis of the fragments; the number of fragments and their masses for each case of metallicity versus rotational energy. 



ivieiaiiiciiy 


Rot. energy 


No. of fragments 


Average fragment mass 


Individual fragment masses 


7/7 


Po 








I0 u 


10 1 


8 


2810.5 


746.9, 1658.2, 2221.0, 3067.8, 3228.8, 3503.6, 3901.4, 4149.8 




10 


5 


5354.6 


1818.4, 5505.1, 6238.3, 6363.6, 6847.7 




10 J 


3 


9363.1 


6047.6, 7196.2, 14845.5 




10" 4 


2 


14075.5 


1 o f r\r\ 1 a a c i o 

13699.2, 14451.8 







2 


14218.3 


7390.8, 21045.8 


10 1 


1 r\-\ 

10 1 


7 


2619.4 


17.8, 1754.1, 2918.7, 3082.4, 3474.2, 3517.1, 3566.8 




1U 


4 


0UjU.4 


3oz4.j, 4/oz.b>, jy/l./, 9jod.d 




10 3 


1 


27235.0 


No fragmentation 




10~ 4 


1 


Ann da r\ 

27849.9 


No fragmentation 







1 


27693.8 


No fragmentation 


10 1 


10 1 




783.0 


53.7, 134.0, 231.0, 256.7, 273.0, 598.4, 632.8, 1374.9, 3492.5 




m-2 
1U 






No fragmentation 




lO" 3 




27067.5 


No fragmentation 




io- 4 




27795.2 


No fragmentation 









28178.1 


No fragmentation 


10" 3 


IO" 1 


9* 


551.0 


28.7, 167.2, 236.1, 308.2, 585.2, 711.1, 821.4, 1037.1, 1064.0 




IO" 2 




24624.6 


No fragmentation 




10" 3 




27122.4 


No fragmentation 




io- 4 




27821.3 


No fragmentation 









28225.5 


No fragmentation 



The two cases marked by a star (*) in the table, represent the upper limits on the number of fragments. 



metallicity then causes y (equation [5]) to be softened to unity or 
less since cooling is more efficient, 



r 



l + 



dlogiT) 
dlogip) ' 



(5) 



In a few low metallicity cases, not discussed in this paper, 
we have not found fragments with densities above the threshold 
criterion. In the same way, the inability to fragment holds for 
all low metallicity runs. While turbulence initially causes den- 
sity perturbations and structures with marginally high density, 
the cloud lacks any development into more compact objects. We 
see that these undeveloped structures are rather diffuse. The in- 
ternal pressure is too high for them to collapse on their own and 
the cooling is not efficient enough to bring down their temper- 
atures. The latter effect causes the thermal Jeans mass, given in 
equation to become very high for these cases and prevents 
further collapse, 



Mj = \ 



5k 



\n) \G\im h ) p 



(6) 



In solar mass u nits, the Jeans equation can be formulated as 
(lFrieswiikll2008h : 



Mj - 7.5( 



;)( 



10K 10 4 cm" 



;T~ 2 We). 



(7) 



Typical Jeans masses are of the order of -10 solar masses 
for a 10-20 K object with a number density of ~10 5 cm -3 . For 
low metallicities, where the temperature can rise up to a few 100 
K, the Jeans mass lies in the range 10 2 - 10 3 M o . Nonetheless, 
the Jeans mass can become as low as a few solar masses at 
higher densities despite the increase in gas temperature. A 
typical Jeans mass diagram of a solar metallicity simulation is 
shown in figure [8] 




-1.0 0,0 1,0 2,0 3.0 4.0 5.0 6.0 7.0 8.0 



log n (cm 3 ) 

Fig. 8. A diagram that shows the Jeans mass as a function of 
number density for a non-rotating solar metallicity run at t = 



From a rotational perspective, we see that faster rotation re- 
sults in more fragments in general. We find fragmentation hap- 
pening for all given initial rotational energies for the highest 
metallicity case, due to the imposed turbulence. Even at the low- 
est rotational energy, collapse to a binary still occurs. Since the 
turbulence imparts significant angular momentum, the j3 = 
case looks similar to the low f3 cases, except that there is more 
diffuse material surrounding the cores. All simulations are su- 
personically turbulent over the first 1-2 dynamical timescales 
and level of to 0.4-2 km/s after 3 ^yn- Typically, the effective 
line widths (Av) at \2t& yn are of the order of 1 km/s on the 
scale of the fragment (~0.1 pc). Higher metallicity tends to lead 
to more turbulent motions. For the lower metallicity cases we 
do not see fragmentation below £c>o%n < 0.03. Here, ojq is 
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Fig. 4. Temperature plots for four different metallicities at t = I2ta yn , Z = Z© (top left), 10 _1 Z o (top right), 10" 2 Z o (bottom left), 
10" 3 Z o (bottom right), for fio = 10 -1 . The color red represents temperatures above 50 K. Contour plots of number density are 
overplotted in these figures in black. Contour levels are 10 6 , 10 5 , 10 4 , 10 3 , 10 2 , 10 1 cm -3 . The axes are in units of 10 18 cm in both 
directions. In the bottom left corner of each image, the Jeans length is illustrated by a white circle. To calculate the Jeans length, 
the average value of the temperatures and densities of the clumped regions of each image is used. These are: T=(10.0, 13.1 13.4, 
19.0)[K] and p=(7.0, 3.3, 0.41, 0.22)[xl0" 19 g cm" 3 ]. 



the angular velocity that is obtained by just putting v = cj^r 
in equation |2l Combining wi th equation [U we get equations [8] 
(Mat sumoto & Hanawall2003l) and[9j In fact, in our studies frag- 
mentation only occurs when the criterion of equation [TO] is met. 
This is found by fitting a function through the points where we 
still see fragmentation; 



foGM 



(8) 



Figure [7J right, shows that the average fragment mass 
<Mf> is higher at lower initial rotational energies. We see here 
that the cloud is fragmenting and the average fragment mass is 
greatly reduced for lower rotational energies when the cloud has 
higher metallicity and is turbulent. 

The specific angular momenta, j, of the fragments with size 
scales of 0.03-0.1 pc are typically a few times 10" 21 cm 2 /s. For 
larger fragments (-0.3 pc) j is so mewhat higher and of the order 
of ~10 22 cm 2 /s in agreement with iBodenheimerl (1 19951) . 



coot dyn exp 



(9) 3.3. The Effect of the Environment on Fragmentation 



' -(log(Z/Z )-1.5) 2 
1.74 



VZ/Z « 1.28 x 10" 



(10) 



We present the density plots of the final stage of every run 
in figure [9] The results of these follow-up simulations can be 
found in tableland the representative average fragment mass is 
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log n (cm °) 



log n (cm °) 




log n (cm 3 ) 



log n (cm 6 ) 





0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 



log n (cm 6 ) 



log n (cm 3 ) 



Fig. 5. The phase diagrams of each metallicity after twelve dynamical timescales for two /3o ratios, 10 1 (left) and 10 3 (right). Top 
to bottom, highest (Z = Z©) to lowest (Z = 10" 3 Z o ) metallicity. 



plotted in figure [TOl 

Dust and CRs in starbursts: When we add an extra heating 
source due to dust, i.e., gas-grain collisional heating, and one 
due to elevated cosmic ray ionization, we see noticeable changes 
in the end results. Initially, the cooling becomes more efficient, 
because the gas has a higher temperature than the dust so that 
the dust acts as a coolant for the gas. When the gas temperature 
drops and falls below the dust temperature, the gas is heated by 



the dust (and cosmic rays). In the end, the fragment temperatures 
are close to the dust temperature of 39 K. The phase diagram for 
a non-rotating case in figure [6] (upper left) shows this steep de- 
crease in temperature, slightly overshooting below the dust tem- 
perature before it settles around the dust temperature. The strong 
collisional coupling between the dust and the gas imposes this 
tight T - n relation. 

Heating by cosmic rays and dust is enough to prevent fur- 
ther fragmentation to smaller scales and causes the largest frag- 
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H 2.0 




0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 



log n (cm 6 ) 



H 2.0 

bJD 



log n (cm 3 ) 



H 2.0 




0.0 i.O 2.0 3-0 4.0 5.0 6.0 7.0 S.O 

log n (cm -3 ) 



H 2.0 



0.0 2.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 



log n (cm 6 ) 



Fig. 6. The phase diagrams of several non-rotating simulations with additional physics for different environments after twelve 
dynamical timescales. The two figures on the right represent isothermal simulations at 25 K (top) and 10 K (bottom). The top left 
one is a phase diagram of the solar metallicity simulation with dust and CR heating. The bottom left figure is for a non-rotating 3D 
simulation also with Z=Z . 




logZ(Z ) log/3 

Fig. 7. Average fragment mass (<Mf>) with initial conditions as determined after twelve dynamical times. Left: <Mf> versus 
metallicity for each initial rotational energy. Right: <Mf> versus rotational energy for every metallicity. 



ments to accrete nearby material. Typically, fewer fragments are 
formed at this point, for the models with fio > 10" 2 . There are 
more fragments for the case fio = 10" 4 where overall gas heat- 
ing due to adiabatic compression and CRs renders T g > Tj. 
This scenario is comparabl e to the situation for population II. 5 
stars dSchneider et al.l2QQ6b . So, dust can enhance fragmentation 
when the conditions are right. 

A remarkable result that we see is that when there is zero 
initial rotation, the warm dusty cloud stops fragmenting in 
contrast to the expected trend. The would-be fragments are 
blurred out and settle in a disk around the sole compact object 
instead. We find that the cloud is very sensitive to initial rotation 
in this regime. A rapid turnover from a fragmenting cloud 



to a non-fragmenting case is found to be around /3o ~ 10 . 
Interestingly, this is in very good ag reement with the lowe r limit 
of fragmentation (J3 = 10 -4 ) that iMachida et alJ (I20Q9I) find 
Starb urst systems may have super-solar metallicities (llsraell 
120091) . We have verified that increasing the metallicity of the 
cloud to twice solar has little impact on our results. 

Isothermal: When we keep the temperature constant, we ef- 
fectively assume that the pressure scales linearly with the den- 
sity. Two near-isothermal runs are compared to the thermal bal- 
ance runs at 10 K and 25 K. In environments that are less ac- 
tive, molecular clouds are known to have these low tempera- 
tures. Although we initialize with these temperatures, accretion 
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Fig. 9. Density plots of the environmental simulations at t = 12^ yn . Molecular clouds in environments having dust and cosmic rays 
are simulated with solar metallicity. Orange to red depicts densities typically higher than 10 4 cm -3 and up to 10 9 cm -3 . Light blue 
portrays the densities smaller than 1 cm -3 . All images have boxsizes of 4x4 parsec. 



shocks, that can increase the temperature up to a few hundred de- 
grees in their wake, still affect the temperature of the cloud and 
the cores. Typically, for gas with density higher than n = 1 cm -3 , 
the temperature remains constant around the initialized value, 
see the two diagrams on the right in figure [6l while the tem- 
perature for the diffuse gas can go up one order of magnitude. 
For the 25 K condition, this makes the Jeans mass higher com- 
pared to simulations with incorporation of the thermal balance 
and where the cooling is dominant. Compared to the highest 
metallicity case Z = Z of the metallicity study, for J3o > 10" 3 , 
the fragments are less well developed and have densities simi- 
lar to the lower metallicity results. Even for the 10 K isothermal 
case, though it has the same final temperature as the cooling in- 
cluded runs, there are fewer fragments. 

At lower rotational energies the difference is less apparent, 
since compressional heating for non-isothermal runs can over- 
come the cooling. In fact, for non-rotating clouds, fragmentation 
for the 10 K isothermal case tends to be higher. This change 
in fragmentation is caused by the fact that there is a tempera- 
ture gradient with increasing density. Although the thermal Jeans 
mass of the metallicity runs is initially higher than in all the near- 
isothermal models, fragmentation is stimulated due to softening 
of y. We note that for the lowest rotational energy (fio = 10" 4 ) 
and the non-rotating case in the solar metallicity model, y is stiff 
(y = 1.4), and the temperature stays relatively high (>100 K) 
compared to the near-isothermal case, because compressional 
heating suppresses fragmentation. The evident difference be- 



tween the two near-isothermal cases (10 and 25 K) can be ex- 
plained by the difference in Jeans mass of a factor ~4. 



Isothermal 25K 
Isothermal 10K 
Dust and CR 




Fig. 10. Average fragment mass versus initial rotational energy 
for simulations within three distinct physical environments at t = 

12fdyn. 
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Simulation 


Rot. energy 


No. of fragments 


Average mass 


Individual fragment masses 




Po 


1 V f 






Dust and CR 


1 r\ 1 

10 1 


1 


1873.9 


454.4, 605.4, 1484.2, 1949.9, 2594.1, 2882.5, 3147.1 




10 


4 


6937.2 


5732.3, 7075.8, 7417.4, 7523.5 




10 J 


3 


9282.3 


6192.2, 6280.7, 15374.0 




10~ 4 


4 


7022.4 


5707.4, 5727.1, 8286.7, 8368.3 







1 


11923.5 


No fragmentation 


Isothermal 25 K 


1 r\-\ 

10 1 


8 


2294.7 


1407.6, 1696.0, 1981.3, 2096.3, 2301.5, 2317.8, 3067.5, 4262.3 




m-2 


J 


O^Jj. / 


ozoo.z, oj/l.J, 11^0/. 4 




io- 3 


3 


9273.5 


7537.8, 8878.8, 11404.0 




IO" 4 


2 


13565.8 


11972.8, 15158.8 







2 


14754.1 


9441.6, 20066.6 


Isothermal 10K 


IO" 1 


7 


3511.7 


1801.5, 1886.3, 2041.8, 3881.1, 4372.6, 5126.8, 5471.8 




IO" 2 


4 


6815.6 


4149.8, 6245.3, 6360.3, 10506.8 




IO" 3 


3 


9302.2 


8452.4, 8939.4, 10514.8 




IO" 4 


3 


9360.1 


5633.7, 11220.1, 11226.5 







3 


16381.3 


9631.5, 23131.2 



4. Conclusions and Discussion 

We have performed 38 hydrodynamical simulations where we 
have rigorously tested the dependence of molecular cloud frag- 
mentation on the ambient conditions as they pertain to starburst 
and dwarf galaxy regions. We have compared these with the well 
known conditions of the Milky Way. Our results on the metallic- 
ity simulations show that the amount of fragmentation and the 
compressibility of the gas scales with increasing metallicity. In 
this, the cooling by gas phase and grain surface chemical species 
plays a subtle but fundamental role. We find that turbulence and 
rotation also play an important role. Turbulence of 0.3£g rav , fol- 
lowing the Larson laws, is a crucial ingredient in the develop- 
ment of fragments. Although rotation is observed for sub-pc 
scale molecular clouds, it is unlikely that clouds of these sizes 
(parsec) have rotational support. All simulations have therefore 
been performed for fio values of much less than unity. These re- 
sults agree with exp ectations and are comp arable to the work of 
iMachidal (12008b and lMachida et alJ (l2009h . albeit at lower den- 
sity. The metallicity has a greater impact on cloud fragmentation 
than Po. For the highest metallicity model, fragmentation (into 
a binary) is still possible for even the lowest (or zero) rotational 
energy. 

Regarding the initial temperature of the simulations, we 
note that starting with a high temperature of 100 K, as for star- 
bursts, does not seem to be too important. When we do not con- 
sider the extra heating sources, radiative cooling can quickly 
overcome the high initial temperature and cool the cloud down 
to -10 K. This is true for most cases except when the initial rota- 
tional energy is very low. In this situation, the adiabatic heating 
and the heating due to shocks, which are more prominent now, 
can overcome the cooling. Fragments have higher temperatures 
>100 K under these circumstances, which suppresses fragmen- 
tation. The surrounding gas, though, stays cold (-10-20 K). 

When we use low metallicity, similar to dwarf galaxies or 
the early universe, the average fragment mass increases in our 
simulations. This leads to massive cloud cores and eventually 
can induce mass ive star formation (Dabring hausen et all 12009: 
lLada et al.|[20Q8h . relevant to the origin of the IMF. 

We summarize our main conclusions as follows: 



- Increased metallicity promotes fragmentation. There is a 
clear and strong dependence of fragmentation on metallicity, 
where the average fragment mass decreases with increasing 
metallicity. 

- Heating by cosmic rays and dust, for high metallicity ini- 
tial conditions, are strong enough to slow down cooling and 
reduce fragmentation. When the gas and dust are thermally 
coupled, we find that the gas temperature stabilizes close to 
the dust temperature, T g = Tj, at densities >10 5 cm -3 . 

- Fragmentation is suppressed for isothermal clouds. In al- 
most all of the near-isothermal simulations that we did, we 
see a smaller number of fragments compared to their non- 
isothermal counterparts, except when adiabatic heating is 
very strong. Although the temperatures are exactly the same 
in several cases, even lower compared to the non-isothermal 
runs for the smallest fio ratios, and the Jeans masses are com- 
parable, the evolution and fragmentation of the cloud differs. 
We attribute this to the stiffer effective equ ation of state, 
i.e. y = 1, that evidently has a great impact dKlessen et al.l 
I2005bl) . 

For the main part in this research, we have performed 2D 
simulations. One can argue that 2D simulations are not represen- 
tative for a 3D cloud. Certain physics, like MHD, is indeed not 
possible to do in 2D, but we have not used any kind of physics 
or initial condition that requires a third dimension. We find that a 
higher merger rate occurs for 2D simulations due to one less de- 
gree of freedom and that our fragments are of somewhat higher 
mass. Still, we do not see many mergers in general. The 3D sim- 
ulations that we have run, show that a disk forms in a time that is 
comparable to or is less than the fragmentation time scale if even 
a modest rotation exists, effectively making the problem a two- 
dimensional one. The results are slightly different, perhaps due 
to the lower resolution, or possibly owing to the shorter simula- 
tion time, but show comparable structures. A disk does not form 
when there is complete lack of initial rotation. This, however, 
does not change our results, since the non-rotating 3D simula- 
tion also fragments into two compact objects with comparable 
masses and enjoys a similar temperature evolution (figure [6]) as 
its 2D counterpart. 

Changes in the density profile do not affect the evolution of 
the cloud as far as the impact of metallicity is concerned. Testing 
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a power-la w density profile (eq.fTTT), as observed in active galac- 
tic regions dGenzel et al.l2003UMapelri et al.l2QQ8l) . against a flat 
profile and keeping the mean density constant, we find that the 
results are similar. The results (not shown) are the same for the 
number of fragments, but also alike in their temperature evo- 
lution. Only the fragment masses are marginall y smaller whe n 
using a flat profile, due to competitive accretion (IBonnelllfeOOSl) . 

po is the central density of n = 10 4 cm -3 , R c is the characteristic 
length of 1 pc. This is the radius within which the central density 
remains approximately constant. R is the specific cloud radius 
and a = 1.4 is the power of the profile. 

We do not reach star densities, but the fragmentation and 
formation of dense cores is determined at the much earlier 
phases that we probe. We have not made use of sink particles, 
as this is usually done to make a jump in resolution and form 
protostars from dense cores, but our fragments should be able 
to form many stars. Opacity effects beyond 10 22 5 cm -2 are not 
considered in our simulations simply because we do not reach 
very high densities. The maximum densities that we reach are 
n~10 8 -10 9 cm -3 . Molecular opacity at high densities can actu- 
ally have a quite different effect on fragmentation compare d to 
densities below 10 3 cm -3 dPoelman & Spaansl l200l 12006b . A 
cloud can become optically thick sooner and trap the lines when 
there are more metals to form molecules with. This in turn will 
heat the system and prevent further fragmentation. Incorporating 
radiative transfer will be the next step in the continuation of this 
research. 
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